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We consider an exclusion process, with particles injected with rate a at the origin 
and removed with rate (3 at the right boundary of a one-dimensional chain of sites. 
The particles are allowed to hop onto unoccupied sites, to the right only. For the 
special case of a = (3 = 1 the model was solved previously by Derrida et al. Here 
we extend the solution to general a, [3. The phase diagram obtained from our 
exact solution differs from the one predicted by the mean field approximation. 
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1 Introduction 



One-dimensional asymmetric exclusion models |l| are of interest for various 
reasons. They are closely related to vertex models @, growth models || 
and, in the continuum limit, the KPZ equation [|J and the noisy Burgers's 
equation. 

Various types of phase transitions occur as a consequence of the interplay 
of particle transport with a localized defect or inhomogenity. Suitably chosen 
boundary conditions can represent the effect of such a defect in an otherwise 
homogenous system. Such transitions have been the focus of many recent 
studies [0,0- @. Some of these models could be solved exactly and allow for 
a detailed study of their steady state properties such as the density profile 
or density correlations |], |7], f| |9j . 

Totally asymmetric simple-exclusion models with nearest neighbour hop- 
ping can be divided into four classes according to the dynamics (sequential 
or parallel) and the choice of boundary conditions (open or periodic). In all 
these models each lattice site i in a chain of N sites is either occupied by 
a single particle (tj = 1) or empty (n = 0) and a particle can hop to the 
neighbouring site in one direction if this site is empty.Q By convention, we 
choose the direction of hopping as to the right. The dynamics can be chosen 
either sequential as in refs. |], |5|, || § or parallel 0, 0, [J. In the case of 
sequential dynamics which we study in this paper particles jump indepen- 
dently and randomly in each time step according to the following rules: At 
each time step t — > t + 1 one chooses at random one pair of sites + 1) 
with 1 < i < N — 1. If there is a particle on site i and site i + 1 is empty, 

then the particle will jump from i to i + 1. All other configurations do not 

1 A model with two different kinds of particles and nearest neighbour hopping has been 
studied in ]iof , a simple-exclusion model where particles can hop over several lattice sites 
in each time step is discussed in [[ll] . 
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change, i.e., 



Ti(t + 1) = Ti(t)r i+l (t) 

(1.1) 

T i+1 {t + l) = T i+1 {t) + {I - T i+l {t))Ti{t) . 

In the case of parallel updating the lattice is divided into neighbouring pairs 
of sites and some stochastic hopping rules are applied in parallel to each pair 
in a first half time step. In the second half time step the pairs are shifted by 
one lattice unit and the same rules are applied again |7|, [|. Both models can 
be defined with either periodic or open boundary conditions. When periodic 
boundary conditions are used, a non-trivial phase diagram can be observed 
by introducing, for example, a single defect || With open boundary 
conditions particles are injected with rate a at the left boundary (which 
we shall call the origin), and absorbed with rate (3 at the right boundary 
05 @) H> HI- Injection and absorption are implemented in the following way: 
when one considers the pair (0,1) where site represents the origin, then the 
occupation number T\{t + 1) of site 1 at time t + 1 is given by 

Tx(t + 1) = 1 with probability ri(t) + a (1 — ri(t)) 

(1.2) 

n (t + 1) = with probability (1 - ot)(l - n(t)) . 

On the other hand, considering the pair (N, N+l) where site N+l represents 
the (right) boundary, the occupation number Tjv(t + 1) at site iV after one 
time step is 

tjv(£+1) = 1 with probability (1 — /3)rAr(t) 

(1.3) 

Tjv(t + 1) = with probability 1 - (1 - /3)r N (t) . 

The model defined by eqs. ( [Lip - ( |1.3| ) can be viewed as a homogeneous 
system connected to a reservoir of fixed particle density a at the origin and 
fixed density 1 — (3 at the boundary. 
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The model with parallel dynamics and periodic boundary conditions with 
a defect was solved by the Bethe ansatz. Recently some steady state 
properties of the model with parallel dynamics and open boundary conditions 
were also found [Dfl. The Bethe ansatz was used also to solve the model with 



sequential updating and translationally invariant periodic boundary 

conditions without defect |12| . The case of sequential dynamics with open 
boundary conditions was studied by Krug || and by Derrida et al 0. 

Krug studied numerically the steady state behavior of this model on the 
line f3 = 1 (Fig. 1). He found, at a = 1/2, a phase transition and an 
associated diverging length scale ||. For a < 1/2 he found an exponential 
decay of the profile to its bulk value with increasing distance r from the 
boundary, while for a > 1/2 the profile decayed as r~ 1 * 2 . Derrida et al [l|] 
expressed the exact steady state and the steady state density distribution for 
arbitrary a and (3 in terms of recursion relations (eqs. (|2.9| ) - (|2.11|) below). 
These recursions were solved explicitly only for a — (3 — 1. For this case 
they showed that the density profile approaches algebraically its bulk value 
Pbuik = 1/2 with increasing distance x from the origin, p — 1/2 ~ x^ 1 ^ 2 . The 
same behavior characterizes the approach of pbuik from below as r -1 ' 2 , with 
increasing distance r from the boundary, confirming the numerical result of 
Krug. 

The phase diagram in the whole a — (3 plane was obtained in |l[ by a mean 
field calculation. Three phases were identified (Fig. 1). In a low density phase 
A, found for a < [3 and a < 1/2, the density profile approaches pbuik = ct 
exponentially with r. This supports Krug's observation of an exponential 
behaviour for (3 = 1, a < 1/2. A high density phase B was found for 
a > (3, (3 < 1/2, which is related to the low density phase by a particle-hole 
symmetry. Here the profile approaches pbuik = 1 — (3 exponentially with x, 
the distance from the origin. Finally, for a, (3 > 1/2 the system is in the 
maximal current phase C. In this phase mean field predicts a power law for 
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Figure 1: Mean field phase diagram of the model in the a — (3 plane as obtained 
in [||. Region A is the low density phase, region B the high density phase and 
region C is the maximal current phase. The phases are separated by the curves 
a = (3 < 1/2 and a = 1/2, /3 > 1/2 and j3 = 1/2, a > 1/2 respectively. The sign 
of the slope of the density profile (as shown in the insets) changes when the line 
a = 1 — (3 is crossed. Note that this is not a phase transition line. 



4 



the profile, with exponent k — 1, whereas the exact result JI| yields k — 1/2. 

Here we present the exact solution to the recursion relations giving the 
steady state and the density profile for arbitrary values for a and j3. We 
show that the phase diagram has a richer structure than that predicted by 
mean field. In particular, we show that there is a phase transition with an 
associated diverging length scale along the two lines, a = 1/2 and (3 = 1/2, 
dividing both the low density phase and the high density phase found in the 
mean field calculation into two different phases. In the low density phase Ai 
defined by a < (3 < 1/2 the profile is exponential, as predicted by the mean 
field calculation. The situation is different, however, in the low density phase 
An , defined by a < 1/2 and f3 > 1/2; there the profile approaches pbuik = ot 
as r _3//2 exp (— r/£) for r ^> 1. This was expected neither from the mean field 
approach nor from the numerical results of Krug. 

The paper is organized as follows. In sec. 2 we present the recursion 
relations obtained in [|] and give an exact solution for arbitrary a and (3. In 
sec. 3 we draw some conclusions from these results and derive the exact phase 
diagram. Then we study the density profile in the various phases for large 
systems (sec. 4) and in sec. 5 discuss in detail the various phase transitions 
that were identified. 

2 Exact solution of the recursion relations 

The steady state of the model defined in eqs. (1) - (3) is given in terms 
of the quantities Pjv( r i> r 2? • • • > t n) which are the probabilities of finding the 
specific configuration of particles represented by the occupation numbers 
(ti,T2, . . • , tn) on the chain with N sites. It turns out to be more conve- 
nient to work with unnormalized probabilities /at(ti, T2, • • • , Tjv) related to 
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Pn(ti,t 2 , . . . ,tn) by 



Pn{ti,t 2 ,...,t n ) = f N (r 1 ,T 2 , . . . ,t n )/Z n (2.1) 

where 

Z n = 12 ••• 12 fN(n,r 2} . . .,t n ) . (2.2) 

71=0,1 T N =0,l 

As shown in [[!]] all the /at(ti,t 2 , . . . , rjy) can be obtained recursively from 
the corresponding quantities in a system with N — 1 sites (eqs. (8) and (9) 
of ref. Q). Here we are interested only in the average occupation number 
(tj) tv of site i in a system of length N, given by 

fa) = T N JZ N (2.3) 

with 

T7V,i = 12 ■■■ 12 T ifN(Tl, T 2 ,..., T N ) . (2.4) 

ri=0,l tjv=0,1 

The normalization and the unnormalized particle density Tjv,i can be 
computed from the quantities 



Yn,k= 12 ■■■ 12 (1-tat)(1 -tat-i) • • • (1 -tx)/jv(ti,t 2 , ...,Tftr) (2.5) 

71=0,1 tjv=0,1 



and 



X V,if = Z • • • Z! (! - T Ar) • • • (1 - T^TpfNin, T 2 ,..., T N ) (2.6) 

71=0,1 TJV=0,1 

by defining 

^v,Af+i = Zn (2.7) 

and 

^"at.jv+i = Tn, p ■ (2.8) 
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(Note that we made a slight change in notation as compared to ref. [JTJ] . 
There the quantities Yn,k were denoted Yn(K) and the quantities Xfj K 
were denoted Xn(K,p).) 

The reason for the introduction of Y^,k for 1 < K < N + 1 and X 7 ^ K for 
£>+l<i^<iV + lis that they can be obtained from the following closed 
recursions |L| 

Y NA = (3Y N _ 1A 

Y N , K = Y NtK ^ + a(3Y N _ ltK for2<K<N (2.9) 

Yn,N+1 = Yn^n + «^Af-l,Af 

with the initial condition 

n,i = p 

(2.10) 

r 12 = a + p . 

These recursions can be simplified somewhat by extending the range of def- 
inition of K to 1 < K < N + 2. If we set K = N + 1 in the second of 
eq. ( |2.9|) , the resulting equation is precisely the third of (|2.9| ), provided we 
use the extended definition Y N _ 1)N+ i = [3~ l Y N _i^ N for [3 ^ 0. Similarly, eqs. 
( p.lOj ) become a consequence of ( [2.9|) by redefining the initial condition as 
Yo,i = 1- These extensions of the definitions of the quantities Y^ k are useful 
in some of the calculations presented below. 

Once the Y N ^ K are determined, the X P N K can be obtained from the recur- 
sion relations ]l[ 

Xnk = X p NK _ x + apX p N _ 1K iorp + 2<K<N 

(2.11) 

X P n,n + i = X p N>N + aX p N _ 1>N forl<p<iV-l 
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with the initial condition 



X N, P +i = aPY N _ ljP+1 for 1 < p < N 



(2.12) 



where we used the extended definitions F ,i an d Yn,n+2 of the Y N ^ K . Solving 
these recursion relations gives the exact average occupation numbers (r.;) 
through eqs. ( |2.3|) , (|2.7|) and (|2.8|) . This was done in ref. |3J for a = [3 = 1. 
Here we present the solution for arbitrary a and f3. 

For a solution of these recursion relations and initial conditions define the 
functions G^ K (x) by 



Af-l 

-<M r„\ _ \ " 



g n,M = X w 



(N>1) 



(2.13) 



r=0 



with 



For later convenience also define £>o,i(0) = = G^ 1 = 1. As a result of the 
symmetries of the coefficients &jv 5 ^(r) these functions satisfy various relations 
given in the appendix. In particular, from the recursion relations ( |A.3| ) and 
the special values ( |A.4[) one can show that the quantity 



K -2 + r 
N 



K-2 



Y N , K (a,(3) = (3 N G N NtK {a) + £ a N ~ s (5 N ~ K+1+s G s ^ N {a) 



(2.15) 



s=0 



solves the recursion relations eq. ( [2.9|) with the initial conditions (|2.10|) . 
Relations (|2.11 ) with initial condition ( 2.12 ) are satisfied by 



N-p 



X p N , K (a,f3) = J2 b N-P,K- P (r)a r+1 (3 r+1 Y N _ r ^ p+1 (a,/3) + 



r=0 



(2.16) 



s 



Using the first of equations ( A.4J) , one obtains from this 



N 

Z N = Y N>N+1 = ]T P'aF-'G'jtbia) (2.17) 

s=0 

and after some computation, involving relabeling of indices, we get 

N-p 

T N , P = X p NjN+1 =af3^2 a s G S N± P ,N-p(P)YN- s -i, P +i(a, (3) . (2.18) 

This expression is exact for any N > 1, 1 < p < N. Substitution in (2.3) 
gives the exact density profile for arbitrary a and (3. 

Eqs. ( |2.17| ) and ( |2.18| ) provide also an exact expression for the conserved 
current j = (r^) — (TjTj + i) = const., and consequently, for the correlation 
function (rjT i+1 ). To see this, note that taking i = N one obtains (t n t n+ i) = 
(1 — (3)(tn) since site N + 1 represents the reservoir of constant density 1 — (3. 
Therefore one has [[L]] 

j = P(t n ) ■ (2.19) 

On the other hand, taking i = one gets (t ti) = a(ri) since site is the 
reservoir of constant density a. Thus we also have Jl| 

j = «(1 - (ti» • (2.20) 

Since, however, our exact result yields (rj) for any i, we can calculate the 
exact current j, and hence (rjTj+i). 

3 Discussion of the density profile 

In order to analyse the density, it is convenient to study the quantity 

tN(p) = (Tn, p +i - T n ^ p )/Z n (3.1) 
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which becomes the spatial derivative of the density profile in the continuum 
limit. It turns out to be given by (p 7^ N) 



t N (p) = (1 - a - (3)P?G%Zl >N - p ((3)a N ^Gl p (a)/Z N . (3.2) 
which for a 7^ 1 — (3 can be more conveniently written in the form 

t N (p) = F p {a)F N _ p {(3)/Z N (3.3) 



with 



and 



F N (x) = x- N - 1 G% N (x) (3.4) 



(1 - a - (3)a N+1 (3 N + 1 



F N {(3)-F N {a) 
a(l - a) - /3(1 - (3) 

1-2/3 



a ^ (3,1- (3 
a = (3^± 



(3.5) 

where the prime denotes the derivative w.r.t. (3. The last of the two eq. ( |3.5| ) 
can be obtained by changing the order of summation in ( p,17[ ). 

For a — 1 — (3 one obtains directly from (|3]2|) that tx(p) = 0, i.e., the 
profile is constant on this curve. This result was already obtained in 
From (|3.3| ) we learn that up to the amplitude Z, the derivative t^ip) of the 
density profile can be written as a product of two functions; one of a and 
the other of (3: t^{p) oc F p {a)F^_ p {f3). 

This fact has important and surprising consequences. It clearly implies 
that phase transitions (i.e. non-analytic changes in the p-dependence of the 
density profile) can occur on two kinds of lines: a = a c and any (3, or (3 = [3 C 
and any a. Hence if a phase transition is predicted to occur on the /3 > 1/2 
segment of the line a = 1/2 (the mean field transition to the maximal current 
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phase), then the transition must extend to the (3 < 1/2 regime as well! This 
means that instead of a single high density phase B, predicted by mean field, 
there are, in fact, two such phases. Indeed, analysis of the function F p (x), 
presented below, reveals that its dependence onp changes at x — 1/2. Similar 
considerations hold for the line f3 = 1/2 and a < 1/2, which separates the 
low density phase A into two distinct phases (Fig. 2) . These new transitions 
were not found by the mean field calculation. 

Another unexpected consequence of the separability into a product is 
the existence of two independent length scales in the model, one determined 
by the injection rate a, the other one by the absorption rate (3. This is 
surprising, as one might believe that only the larger of these two quantities 
determines the behaviour of the system. In fact, as long as the system is not 
in the maximal current phase, this indeed is the far as the current 

j = (rj) — (rjTj+i) is concerned: In the continuum limit one has j = (3(1 — j3) 
for a > j3, [3 < 1/2, and j = a(l — at) for (3 > a, a < 1/2, (whereas 
j = 1/4 if both a and (3 are larger than 1/2). Since in the mean field 
calculation the shape of the density profile is determined by only the current, 
phase transitions are seen neither at a = 1/2, (3 < 1/2, nor at (3 = 1/2, 
a < 1/2. Prior to presenting an explanation for the unexpected existence of 
the additional phases and phase transitions, we study the density profile in 
the thermodynamic limit iV — >• oo. 

4 Density profile in the large N limit 

We want to discuss the density profile of a large system (N ^> 1) as a function 
of the space coordinate p, at large distances from both ends, i.e., we consider 
p ^> 1 and r = N—p ^> 1. So we need an asymptotic expression for Fl(x) for 
large L. Splitting Fl(x) into two pieces F^\x) and F < f'\x) as in eq. (|A.11|) 
allows for an expansion in 1/L. For x < 1/2 the dominating contribution is 
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Figure 2: Exact phase diagram of the model in the a — (3 plane. The low (high) 
density phase shown in Fig. 1 is divided into two phases Ai and An (Bi and Bn) 
along the curve (3 = 1/2 (a = 1/2). 
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since /F^ oc exp (—aL) with some constant a: 
1 - 2x 



Fl(x) = ; 1 + O(e^) x<\. (4.1) 

If x > 1/2 then F^ = and up to order 1/L 



4 L 



(4.2) 



X > i 



~ (l-2x) 2 v^FL3/2 + — 2 

This expression diverges for x — > 1/2. For x = 1/2 one obtains (see ( |A.14j )) 



F L (x) = 2^= l+OOL- 1 ) . (4.3) 



Using the expansions eqs. Q4.1|) - ( |4.3|) and the expression fl3.5|) for the 



normalization Z^r one can compute the shape of the density profile given by 
t^{p). We define a length scale £ CT by 

C^-logCMl-a)). (4.4) 

As cr reaches 1/2, £ CT diverges. For the various phases Aj - C (Fig. 2) one 
finds in the large N limit (such that 1 <C p, 1 <C iV - p, i.e., p is far from 
both ends of the system) the following results: 

High density phase E>i : 

This phase is defined by the region (3 < a < 1/2. From the expansion (JO) 
one finds an exponential decay of the density profile with a exponential decay 
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with length scale £ 1 = £ a 1 — £ 1 



/9 ' 



™-<'-->KH)(£H)' 

(4.5) 

= (1 -2q) W 
v ' V 4a(l - a) / 



The density approaches its bulk value pbuik = 1 — (3 from below. One has 
j = (3(1 - (3) and from (g^Ol) (n) = 1 — — /9)/a < 1 - (3 = (r N ). 



Transition line from high density phase Bj to high density phase B n 



On approaching a = 1/2 from below in the region (3 < 1/2, we find that £ a 
diverges but £g remains finite. For a = 1/2 the slope of the profile is of the 
form 

t N (p) ~ p-*» e - p/ fr . (4.6) 
The values of the length scale ^ and the exponent be read off the 



exact expression (fOf) which for large N becomes 

(1 -2/?) 2 (4/3(1 -/?))»> 

Mp) " ^ — 

(4.7) 

_ i l ' 2 (3) -1/2 -p/fr 

2^ P 

The current and the boundary values are given by the same expressions as 
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in the high density phase I. 



High density phase Bn : 

On crossing the phase transition line into the high density phase Bn defined 
by a > 1/2 and (3 < 1/2, the decay exponent changes to z a = 3/2 (see eq. 
(fO)) and one obtains 



(1- a -/?)(<*-/?) (4/3(1-/?))* 
tN{P) ~ (l-2a)V5F ^ 

= 3/2 -p/ C „ 

(1-2«) 2 V ^F P 



(4i 



The current and the boundary values are given by the same expressions as 
in the high density phase I, but note that the slope of the profile changes 
sign on the curve a = 1 — (3. Along this curve the density is constant, 
( T i) = Pbuik = 1 — P for 1 < i < N. For a > 1 — (3 the slope is negative. 

Transition from high density phase Bn to the maximal current phase C: 



When (3 reaches the critical value 1/2 in the region a > 1/2, (3 < 1/2, then 
also diverges and the slope of the profile is given by 



Near the origin (1 <^ p <^ N) we can neglect the piece with p/N in ( f4.9| ), 
so the slope is dominated by p~ Za with z a = 3/2. In the boundary region 
(p = N — r, 1 C r < N) the shape of the profile is determined by r~ z < 3 with 
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zp = 1/2, but the amplitude of t N (r) is only of order 1/N. Therefore, up 
to corrections of order 1/N, the profile near the boundary is flat, whereas it 
decays as p' 1 / 2 with the distance p from the origin to its bulk value p bu ik = 
1/2. The The current reaches its maximal value j = 1/4 and one finds 
(tn) = Pbuik = 1/2 and (n) = 1 - l/(4a). 



Maximal current phase: 



If (3 > 1/2 and a > 1/2, the derivative t^(p) depends neither on a nor on 
p. Near the origin and near N the slope of the profile is determined by 
z a = zp = 3/2: 

t N ( P ) = --L(i - P /iv)- 3 V 3/2 (4.io) 

4V 71 " 

Therefore the density approaches its bulk value p bu i k = 1/2 as p' 1 / 2 with 
the distance p from the origin from above and as r~ x l 2 with the distance 
r = N — p from the boundary from below. The current takes its maximal 
value j max = 1/4 throughout the phase and one obtains (t n ) = 1/(4/9) and 
(n> = 1 - l/(4a). 



Low density phase Ai : 



This phase is defined by a < (3 < 1/2 and is related to the high density phase 
Bi by a particle-hole symmetry and therefore the decay is exponential. One 
finds 



t N (p) = (1-2/3) 1- 



4afl - a)\ {4a(l-a)\ N p 



(4.11) 



e-r/C 
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with a length scale £ _1 = — £1 and r = N — p 1. The density 
approaches its bulk value pb u ik = 01 from above. The current is given by 
j = a(l — a) and therefore according to (|2.19|) (tn) = a(l — a)//3 > a = (ti). 



Low density phase An : 

The profile in this regime {(3 > 1/2, a < 1/2) is obtained from (|4.8|) by 
exchanging a and /3 and substituting p by r = N — p. This is a result of the 
particle-hole symmetry of the model. In the same way one obtains the profile 
on the phase transition lines from Aj to An and from An to C out of the 
profiles on the phase transition lines from Bi to Bn and Bn to C respectively. 



Coexistence line: 



If a — (3 < 1/2 both £ a and ^ are finite, but since £ a = one gets £ 1 = 0. 
As a result one finds a linear profile with a positive slope 

t N (p) = (1 - 2a) /N . (4.12) 

The current is given by j = a{l — a) and one has (ti) = a and (rjv) = 1 — a. 

5 Discussion of the phase diagram 

We turn now to discuss the various phases and the transitions between them 
on a more physical, intuitive basis. First we consider the case (3=1. This 
situation corresponds to connecting the system to a reservoir of fixed density 
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Po = ot at the origin, and another "reservoir" with pN+i — 1 — P — at the 
boundary. We will consider the limit iV — > oo, and ask what are the possible 
steady state density profiles that the system can have, and which interpolate 
between the two limiting values p an d Pn+i- 

Let us start with a < 1/2, and try a density profile (a) that approaches 
(for 1 < x < Af) a constant bulk value p < a, before it decays to pn+i = 0. 
We now show that such a profile cannot be a steady state. To see this, note 
that in a bulk region with constant density there are no correlations (the 
steady state factorizes into a product measure) and therefore the current in 
the bulk is given by j = p(l — p); whereas at the origin it is jo = a(l — pi), 
where pi, the density at x = 1, satisfies 1/2 > a > p\ > p. If we can 
show that jo > j, particles accumulate between the origin and the bulk, 
and hence the density is not stationary. Clearly, for pi = p we have jo = 
a{l — p) > p(l — p) = j since a > p. On the other hand, for pi = a we have 
jo = a(l — a) > p(l — p) = j for 1/2 > a > p. Hence j > j at the two 
endpoints of the interval [p, a] to which pi is limited; and since jo is a linear 
function of pi, we must have jo > j for the entire interval. 

A different possible steady state profile (b) is one with a < p < 1/2. 
Here we can show that jo = a(l — a) < p(l — p) = j: Under the present 
assumptions the density first interpolates between a and p, and hence a < 
pi < p, and as before, the relationship we wish to prove holds at both 
endpoints of this interval. If this holds, however, more particles leave the 
bulk than enter it, and p must decrease. Thus also (b) cannot be a steady 
state. 

The last possibility of the kind considered, (c), has p > 1/2 > a; we will 
return to this case later and show that for the presently used values of a and 
f3 it cannot be a steady state profile either. The only remaining situation is 
the one in which p = a. Then, obviously, jo = a(l — a) — j. Hence the bulk 
steady state density must equal that of the reservoir. 
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The assumption a < 1/2 was crucial for our proof of this fact, which is 
no longer true if a > 1/2. In that case the bulk density is, independently of 
a, given by p = 1/2. To see this, we again assume all other possible values 
for the bulk steady statedensity, and rule out every other scenario. Let us 
start by assuming a decay to a bulk value p < 1/2; since, supposedly, we are 
in a steady state, we can choose some point i at which l/2>pj>pasa 
new initial point of fixed density a' = pf, a' plays now the role of a < 1/2 
of the previously discussed situation (a'), which, as we have shown, cannot 
be a steady state. Another possibility (&') has a > p > 1/2. In this case we 
recall that near N the density profile must go from p to pn+i = 0. To rule 
this out, we view the site N + 1 as a reservoir of holes of fixed density 1. The 
bulk with p > 1/2 corresponds to hole density p^ — 1 — p < 1/2; holes move 
to the left, and if we exchange the roles of holes and particles, this situation 
becomes precisely the case (a') discussed above. Hence (b') is not possible 
either. | 

We have just shown that for a > 1/2 no steady state is possible with 
either bulk density p < 1/2 or p > 1/2; hence the only possibility left is 
Pbuik = 1/2. That is, for a > 1/2 the bulk density is that one which supports 
the maximal current, irrespective of a, the density of the reservoir. This 
explains the transition observed at a = 1/2, from a low density phase with 
Pbuik = a to the maximal current phase, in which pi, u ik = 1/2. 

For the sake of convenience we limited the previous discussion to the 
(3 = 1 line. We now show that the transition survives when we move off this 
line. Of the cases discussed above, (a), (b) and (a') were ruled out with no 
mention of the fact that (3 = 1. In case (&') (and its equivalent, case (c)), 
we used a particle-hole symmetry to map the situation onto case (a'). Since 
there we had a > 1/2, the same argument goes through for holes if (3 > 1/2; 

2 Note that the situation (c) of the a < 1/2 case, to which we promised to return, also 
requires that p goes from p > 1/2 to pn+i = 0, and therefore is ruled out in the same way 
as (&')■ 
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hence the same considerations give rise to the same phases as were obtained 
on the (3=1 line. This completes the picture for the regions 1/2 < (3 < 1, 
and by particle-hole symmetry, for l/2<a<las well. 

Note that in the low density phase An with pbuik = a the slope of the 
density profile changes sign on the curve a — 1 — (3. This can be understood 
as follows. The probability that a particles moves in the bulk (its average 
velocity) is v — 1 — a, while the probability that it moves at the boundary 
is vn = P- If /? > 1 — a then v n > v and the system becomes depleted near 
the boundary because the current j = pv is conserved. This corresponds the 
negative slope of the profile in this regime. On the other hand, if (3 < 1 — a 
one has vn < v and particles pile up. This leads to a positive slope. In 
the high density phase Bn one finds the same result when compairing the 
velocity v o = 1 — a at the origin with the bulk velocity v — (3. 

Next, we discuss the low density phase Ai , the high density phase Bi 
in the region < at, (3 < 1/2 and the transition between them. The bulk 
values in both phases and the slope of the profile can be derived in the 
same way as for the phases An and Bn . Note that in both phases one has 
vq = 1 — a > j3 = vn and therefore the slope is always positive (particles pile 
up). The current is given by 

j = min (a(l - a), (3(1 - (3)) . (5.1) 

In order to understand the shape of the profile in phases Ai and Bi we assume 
that it is built up by a superposition of profiles with a constant density a up 
to some point Xq, followed by constant density 1 — (3. We call this sudden 
increase of the average density a domain wall since it separates a region 
of high density 1 — (3 from a region of low density a. The picture we have 

3 We assume the width of this domain wall to be very small compared to the size of the 
system. 
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in mind for this scenario is that particles injected with rate a at the origin 
move with constant average velocity 1 — a > 1/2 until they hit the domain 
wall where they get stuck and continue to move only with velocity (3 < 1/2. 
This region of high density is caused by the blockage introduced through the 
connection to the reservoir of density 1 — (3 at the boundary. Such a scenario 
is plausible, since constant densities a < 1/2 starting from the origin and 
1 — (3 > 1/2 connected to the boundary are both stable situations of the 
system as discussed above. The probability p(x) oc exp (— x/C,) of finding 
this domain wall at position x is determined by the length scale £ given 
by f- 1 = - (see O and (|X[fl). If a < (low density phase) 
then particles are absorbed with a higher probability than they are injected 
and the probability of finding the domain wall decreases exponentially with 
increasing distance r = N — x from the boundary. On the other hand, in 
the case where a > (3 (high density phase) the situation is reversed and p(x) 
decreases with increasing distance from the origin. Averaging over all such 
profiles with the weight p(x) leads to the observed exponential decay to the 
respective bulk value. This picture provides also a natural explanation of the 
linear profile on the transition line a = (3 where the absorption and injection 
probabilities are equal. Here £ diverges and the probability of finding the 
domain wall at x is independent of x. Averaging over step functions with 
an equal weight for every position of the step clearly gives a linear profile. 
It is worth noting that the mean field calculation done in |IJ gives a correct 
description of the phases Ai and Bi , but it singles out the constituent step 
function with the domain wall located in the center as the profile on the 
phase transition line. We should mention that this analysis, in particular 
the use of the domain wall picture for a description the two phases and the 
phase transition line, is based on our studies of a similar exclusion process 
with open boundary conditions but parallel dynamics f9j. For this model 
we found phases of type Ai and Bi and a phase transition separating them 
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as in the system with sequential dynamics studied here. A careful study 
of the equal time correlation functions leads to our interpretation in terms 
of domain walls. As the transition lines to the the phases An or Bn are 
approached, this picture becomes invalid. 

Finally we briefly discuss the phase transition from the high density phase 
Bi to the high density phase Bn . On approaching a = 1/2 (but (3 < 1/2) 
the length scale £ a diverges while £p remains finite. As a result neither the 
bulk density nor the way how the bulk density is approached depends on a 
(except for the trivial fact that the density at the origin and consequently 
the amplitude of the derivative of the profile depend on a). The decay to 
the bulk density = 1 — (3 is determined by £g alone. Similarly the 

current does not depend on a, being j = (3(1 — (3). A description of phase Bn 
also in terms of constituent profiles is appealing at first sight, but it is less 
convincing since a constant profile of density a > 1/2 at the origin is not a 
stable situation. In order to get a more intuitive insight regarding this phase 
transition, we consider again the transition from the low density phase An to 
the maximal current phase on the line a = 1/2 but (3 > 1/2. In the maximal 
current phase C the bulk density and the way how it is approached does not 
depend on a whereas in the low density phase An a does determine the bulk 
density and how the profile decays to it. This is obviously due to the fact 
that if a > 1/2 the particles close to the origin block each other rather than 
flowing away. As a result, the information corresponding to a change in the 
injection rate does not penetrate into the system. Clearly this description of 
the effect of a increasing beyond 1/2 on the transition from the low density 
phase Ai to the maximal current phase does not depend on the absorption 
at the boundary and is therefore also applicable to the transition from phase 
Bi to phase Bn .0 

4 As a result of the particle-hole symmetry of the problem, the discussion of the transi- 
tion from the low density phase Ai to the low density phase An proceeds along analogous 
lines. 
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We conclude that the phase transitions to the maximal current phase 
from the phase An (or Bn ) is of the same nature as the phase transition 
from Bi to Bn (or from Ai to An ) . These transitions are caused by reaching 
the maximal transport capacity of the system at the origin (or boundary) 
and result in the divergence of the corresponding length scale determining the 
shape of the profile. Note that in our explanation it was necessary to take into 
account local correlations rather than only the current. This is the reason 
why these phase transitions are not found in the mean field calculation. 

As opposed to these transitions, the one at a — (3 that takes the system 
from the low density phase Ai to the high density phase Bi is caused by the 
building up of domain walls. Such a wall is generated by the inhomogenity 
forced on the system by being connected to two reservoirs of different densi- 
ties. At the transition line the wall can be anywhere with equal probability. 
The "coexistence" of low and high density regions is known to occur also in 
other systems with such an inhomogenity || [7|, |9|, [T0| . 



Note added: After completion of this work we received a preprint by 
B. Derrida, M.R. Evans, V. Hakim and V. Pasquier, who solved the same 
problem by a different method. 
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Appendix Some useful identities for the G- 

function 



The function Gf} K (x) was defined in eq. (|2.13|) as 



with 



Af-l 

G^ K {x) = £ b N , K (r)x r (N > 1) (A.l) 

r=0 



l / \ ( K-2+r\ ( K-2+r\ , . nS 

WW= K-2 )-[ N ■ ^ 



Furthermore we defined 6 ,i (0) = G\ = Gg 1 = 1. Using the symmetries of 
the coefficients h^ t jc(r) it is easy to prove the following recursion relations: 

G^ K {x) = G%_ ltK (x) + xGfa^x) (2<K<N+1) 
G N N , K {x) = G^ix) + xG^zi K (x) (2<K<N + 1) 

G% N (x) = GfJ^^+xG^ix) (N>2) 



(A.3) 



One also finds 



Gn,n( x ) 


— G N ^ N+1 (x) 


— G^^n ( x ) 


= G%$ +1 {z) (N>1) 


G°n } k( x ) 


= G^ N+2 (x) 







Gn } k( x ) 


= G% tl (x) 


1 


(K, M < N + 1 



(A.4) 

Eqs. (|A.3|) and (|A.4j) are necessary to prove that the function Y^xiot, 0) 
and Xff K (a, (3) satisfy the recursion relations and initial conditions Q2.9| ) - 

&■ 

Gjl k( x ) can be expressed in terms of incomplete /3-functions: 

{l-x) K ~ l G K l K {x) = I 1 _ x (K-l,M)-(^-—j I^ x (N+l,M-N+K-2) 

(A.5) 

(This is a direct consequence of the definition of I X (P, Q).) From this expres- 
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sion one obtains 



(1 x) N+1 G^ K (x) x M ' G^\ M _ N+K _ 1 {\ x) — (1 x) K+2 x N K+2 

(A.6) 

and by setting M = N and M = N + 1 



G^ K 2 (x) = (1 - x)G»XIk + i(x) (2 < K < N + 1) . (A.7) 



From this relation one can see that relations ( |A.3| ) are consistent for M = 
N = K. 

In the expression ( |3.3| ) for the density profile only the function G^ L (x) 
appears. From eq. ( |A.7| ) one finds 

G L L Ax) = (1 - x)G L L f 1 r +1 {x) + c L x L+1 (A.8) 



with 



c, = -^(L + l) = z J^ . (A.9) 
Defining Fl(x) = x^~ L ~ l 'G^ L (x) one gets 

Fl(x) = 1 = * +1 - Ec. (x(l - x)) k ~ L . (A.10) 

— X)) k=0 

Now we can reexpress Fl (x) for x ^ 1 /2 in terms of a hypergeometric func- 
tion: 

(1 - 2x)G(l - 2x) 



(x(l — x)) 



(l-2x) 2 ~ ' "- 2 " J 

F£\x)+FP (x) 
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Here Q(z) is the step function 

' 1 z > 



e(z) = 



(A.12) 

^ < 



For large L the hypergeometric function F(l, |; L + 2; z) reduces to 

F(l, §; L + 2; z) = 1 + . (A.13) 

Special exact values of Fl(x) are 

^) = 2( 2 L L ) = ^^(l + 0(L- 1 )) (A.14) 



and 

4 L 



F L (l) = c L = - m ^(l + 0(L-i)) . (A.15) 
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